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Abstract. A self-consistent electrodynamics of a parti- 
cle accelerator in a rotating neutron-star magnctosphcrc 
is investigated on the two-dimensional poloidal plane. 
Solving the Poisson equation for the electrostatic poten- 
tial together with the Boltzmann equations for electrons, 
positrons and gamma-rays, it is demonstrated that the 
created current density increases to be super-Goldreich- 
Julian if the trans-field thickness of the gap becomes 
thick enough. This new solution exists from the neutron- 
star surface to the outer magnctosphcre with a small- 
amplitude positive acceleration field in the inner part, 
which works to extract ions from the stellar surface as a 
space-chargc-limitcd flow. The acceleration field is highly 
unscreened in the outer magnetosphere, in the same man- 
ner as in traditional outer-gap models. 



1. Introduction 

The Energetic Gamma Ray Experiment Telescope 
(EGRET) aboard the Compton Gamma Ray Observa- 
tory has detected pulsed signals from at least six rotation- 
powered pulsars. Since interpreting 7-rays should be less 
ambiguous compared with reprocessed, non-thermal X- 
rays, the 7-ray pulsations observed from these objects are 
particularly important as a direct signature of basic non- 
thermal processes in pulsar magnctospheres. 

Attempts to model the particle accelerator have tradi- 
tionally concentrated on two scenarios: Polar-cap models 
with emission altitudes within several neutron star radii 
over a polar cap surface (Harding, Tademaru, & Esposito 
1978; Daugherty & Harding 1982, 1996; Stumer, Dermer, 
& Michel 1995), and outer-gap models with acceleration 
occurring in the open zone located near the light cylin- 
der (Cheng, Ho, & Ruderman 1986a, b; Chiang & Romani 
1992, 1994; Romani and Yadigaroglu 1995). Both models 
predict that electrons and positrons are accelerated in a 
charge depletion region, a potential gap, by the electric 
field along the magnetic field lines to radiate high-energy 
7-rays via the curvature and inverse-Compton processes. 

The outer gap models have been successful in explain- 
ing the observed light curves, particularly in reproducing 



the wide separation of the two peaks commonly observed 
from 7-ray pulsars (Kanbach 1999; Thompson 2001). In 
these models, they consider that the gap extends from the 
null surface, where the Goldrcich- Julian charge density 
vanishes, to the light cylinder, beyond which the velocity 
of a co- rotating plasma would exceed the velocity of light, 
adopting the vacuum solution of the Poisson equation for 
the electrostatic potential. 

However, it was analytically demonstrated by Hi- 
rotani, Harding, and Shibata (2003, HHS03) that the gap 
inner boundary shifts towards the star as the created cur- 
rent increases and at last touch the star if the created 
current exceeds the Goldreich- Julian (GJ) value at the 
surface. Therefore, to understand the particle accelerator, 
which extends from the stellar surface to the outer mag- 
netosphere, we have to merge the outer-gap and polar-cap 
models, which have been separately considered so far. 

In traditional polar-cap models, the energetics and pair 
cascade spectrum have had success in reproducing the ob- 
servations. However, the predicted beam size of radiation 
emitted near the stellar surface is too small to produce the 
wide pulse profiles that are observed from high-energy pul- 
sars. Seeking the possibility of a wide hollow cone of high- 
energy radiation due to the flaring of field lines, Arons 
(1983) first examined the particle acceleration at the high 
altitudes along the last open field line. This type of ac- 
celerator, or the slot gap, forms because the accelerating 
electric field, En, is screened at increasingly higher alti- 
tude as the magnetic eolatitude approaches the edge of the 
open field region (Arons & Scharlemann 1979). Muslimov 
and Harding (2003, 2004a,b, hereafter MH03, MH04a,b) 
extended this argument by including two new features: ac- 
celeration due to space-time dragging, and the additional 
decrease of E\\ at the edge of the gap due to the narrow- 
ness of the slot gap. It is noteworthy that the polar-slot 
gap model proposed by MH04a,b is an extension of the 
polar-cap model into the outer magnetosphere, assuming 
that the plasma flowing in the gap consists of only one 
sign of charges, which leads to a negative E\\ when the ro- 
tation and magnetic moment vectors resides in the same 
hemisphere. However, we should notice here that the elec- 
tric current induced by the negative En contradicts with 
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the global current patterns, which is derived by the EMF 
exerted on the spinning neutron-star surface, if the gap is 
located near the last-open field line. 

On these grounds, we are motivated by the need to 
contrive an accelerator model that predicts a consistent 
current direction with the global requirement. To this 
aim, it is straightforward to extend recent outer-gap mod- 
els, which predict opposite Eu to polar-cap models, into 
the inner magnetospherc. Extending the one-dimensional 
analysis along the field lines in several outer-gap models 
(Hirotani and Shibata 1999a, b, c; HHS03), Takata, Shi- 
bata, and Hirotani (2004, hereafter TSH04) and Takata 
et al. (2006, hereafter TSHC06) solved the Poisson equa- 
tion for the electrostatic potential on the two-dimensional 
poloidal plane, and revealed that the gap inner boundary 
is located inside of the null surface owing to the pair cre- 
ation within the gap, assuming that the particle motion 
immediately saturates in the balance between electric and 
radiation-reaction forces. In the present paper, we extend 
TSH04 and TSHC06 by solving the particle energy distri- 
bution explicitly, and by considering a super-GJ current 
solution with ion emission from the neutron star surface. 



2. Basic Equations and Boundary Conditions 

We formulate the basic equations to describe the particle 
accelerator, extending the method first proposed by Be- 
skin ct al. (1992) for black-hole magnetospheres. The first 
kind equation we have to consider is the Poisson equation 
for the non-corotational potential, 
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where the general rclativistic Goldrcich- Julian charge den- 
sity is defined as 
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The second kind equations we have to consider is the 
Boltzmann equations for particles. Imposing a stationary 
condition d t + Sld,p = 0, utilizing V B = 0, and intro- 
ducing dimensionless particle densities per unit magnetic 
flux tube such that n± = N±/{Q,B /2nce), we obtain the 
following form of the particle Boltzmann equations, 



ccosx- 



dn± dp dn± d\ dn± 



ds 



dt dp dt dx 
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where ,s denotes the distance along the magnetic field line, 
(p,X) the momentum and pitch angle of the particles. The 
upper and lower signs correspond to the positrons (with 
charge q = +e) and electrons (q = — e), respectively; p = 
\p\ and ds/dt = ccosx, 



dp „ Psc 

— = qE\\ cos x 

dt " c 



(4) 



d X qE\\smx dilnB 1 /' 2 ) 



dt 



sin*, 



(5) 



For outward- (or inward-) migrating particles, cosx > 
(or cosx < 0). Since we consider relativistic particles, we 
obtain T = p/(m e c). Using n±, we can express p c as 

£113 f f 

Pc= 9^l t n +( s ' 6 '*' r '^) ~ n -( s '°*i T >x)]drdx- (6) 

The synchro-curvature radiation-reaction force, Psc/ C , is 
given by Cheng and Zhang (1996). Collision terms S± rep- 
resents the particle appearing and disappearing rate from 
(p,x) and are given by the inverse-Compton scatterings 
and one-photon and two-photon pair creation processes. 

The third kind equations we have to consider is the 
Boltzmann equation for 7-rays. Imposing the stationary 
condition, we obtain 
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where (p = tp — Ot. To compute k l 1 we have to solve the 
photon propagation in the curved spacetime. For the com- 
putation of k l and 5 7 , sec Hirotani (2006). 

In order to solve the set of partial differential equa- 
tions |T]), ([3]), and (J?]) for n±, and g, we must impose 
appropriate boundary conditions. We assume that the gap 
lower boundary, 8* — (9™ ax , coincides with the last open 
field line. Moreover, we assume that the upper boundary 
coincides with a specific magnetic field line and parameter- 
ize this field line with = #™ m . Determining the upper 
boundary from physical consideration is a subtle issue, 
which is beyond the scope of the present paper. There- 
fore, we treat #™ m as a free parameter. We measure the 
trans-field thickness of the gap with 



h m = 



(8) 



If h m = 1.0, it means that the gap exists along all the open 
field lines. On the other hand, if h m <C 1, the gap becomes 
transversely thin. To describe the trans-field structure, we 
introduce the fractional height as 



(9) 



Thus, the lower and upper boundaries are given by h = 
and h = /i m , respectively. The inner boundary is assumed 
to be located at the neutron star surface. For the outer 
boundary, we solve the Poisson equation to the light cylin- 
der (see Hirotani 2006 for details). 

First, to solve the elliptic- type equation (Q}, we impose 
$ = on the lower, upper, and inner boundaries. At the 
outer boundary, we impose d^J /ds = 0. Generally speak- 
ing, the solved E\\ = —{d^ /ds) s ^o under these boundary 
conditions does not vanish at the stellar surface. For a 
super-GJ current density in the sense that p e — pqj < 



106 



K. Hirotani: Particle Accelerator in Pulsar Magnetospheres 



-i 1 1 r 



i 1 1 r- 




Fig. 1. Created current density j c (in unit of QB/2ir) as 
a function of the transfield thickness h for five different 
trans-field heights. 

holds at the stellar surface, equation (JTJ) gives a positive 
electric field near the star. In this case, we assume that 
ions are emitted from the stellar surface so that the ad- 
ditional positive charge in the thin non-relativistic region 
may bring Eu to zero. 

Secondly, to solve the hyperbolic-type equations ([3]) 
and ([7]), we assume that neither positrons nor 7-rays are 
injected across the inner boundary; thus, we impose 

n+(s in A,r, X )=0, 5 (s in A,£ 7 ,0 7 ) = O (10) 

for arbitrary 9*, T, < \ < tt/2, E 1i and cos(0 7 — 9b) > 0, 
where 9b designates the outward magnetic field direction. 
In the same manner, at the outer boundary, we impose 

n_( s ° ut A,r, x ) = o, 5 ( s ° ut A,£ 7 ,0 7 ) = o (11) 

for arbitrary T, tt/2 < x < 7r, Ey, and cos(£? 7 — 9b) < 0. 

3. Application to the Crab Pulsar 

We apply the scheme to the Crab pulsar. We present the 
created current density, j Cl as a function of h, for the 
five different trans-field thickness, h m = 0.047 (dotted), 
0.048 (solid), 0.060 (dashed), 0.100 (dash-dotted), and 
0.160 (dash-dot-dot-dotted). The thin dashed line repre- 
sents |pgj/(^-B/27tc)| s= o; if jo appears below (or above) 
this line, the created current is sub- (or super-) GJ along 
the field line specified by h. For h m = 0.047, the solution 
(dotted curve) is sub-GJ along all the field lines; thus, 
screening due to the discharge is negligible as shown by 
the dotted curve in figure which depicts Eu(s, h) at the 
central height h = h m /2. As h m increases, the solution be- 
comes super-GJ from the higher latitudes, as indicated by 
the solid, dashed, dash-dotted, and dash-dot-dot-dottcd 
curves. As a result, £j| is screened significantly with in- 
creasing /i m , as demonstrated by figure[H This screening 
of E\\ has a negative feed back effect in the sense that j e 
is regulated below unity. 




Fig. 2. Field-aligned electric field at h = h m /2 as a func- 
tion of s/tolc- The dotted, solid, dashed, and dash-dotted 
curves corresponds to h = 0.047, 0.048, 0.060, and 0.100, 
respectively. 

To understand the screening mechanism, it is helpful 
to examine the Poisson equation ((T|). In the transversely 
thin limit, it gives 



D.B* 2-kc 
c VLB 



>)• (12) 



We thus find that E^ = —d^/ds is approximately propor- 
tional to —d(p/B — pqj/ B)/ds. It is, therefore, important 
to examine the two-dimensional distribution of p/B and 
Pqj/B to understand £j|(s,/i) behavior. 

In figure [3l we present pj (HB /2ttc), p c /(flB /2ttc), 
and pgj/{VLB/2-kc), as the solid, dash-dotted, and dashed 
curves, at nine discrete magnetic latitudes at h = (3/8)/i m , 
h m /2, (5/8)/i m , (3/4)/i m , for the same parameters as fig- 
ure O If there is a cold-field ion emission from the star, 
the total charge density (solid curve) deviates from the 
created charge density (dash-dotted one) . The current be- 
comes super-GJ for h > (6/16)/i m ~ 0.0225. Along the 
field lines with super-GJ current, p c — pqj becomes neg- 
ative close to the star. This inevitably leads to a positive 
Eu , which extracts ions from the stellar surface. 

In the outer region, p/B levels off in s > 0.5wlc for 
h > h m /2. Since pcj/B becomes approximately a linear 
function of s, .Em remains nearly constant in s > 0.5iz7lc, 
in the same manner as in traditional outer-gap model. In 
the inner region, on the other hand, inward-directed 7-rays 
propagate into the convex side due to the field line curva- 
ture, increasing particle density with h. A two-dimensional 
effect in the Poisson equation becomes important in the 
higher altitudes (O.Itolc < s < 0.3h7lc) along the higher- 
latitude field lines(/i > 0.0225). Outside of the null surface, 
s > 0.09tz7lc, there is a negative p e s in the sub-GJ current 
region (h < 0.015). This negative p c g works to prevent En 
from vanishing in the higher latitudes, where pair creation 
is copious. However, the created pairs discharge until E\\ 
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vanishes, resulting in a larger gradient of p than that of 
Pgj hi the intermediate latitudes in 0.0225 < h < 0.0262. 
In the upper half region (0.03 < h < h m = 0.06), dp/ds 
does not have to be greater than dpcj/ds, in order to 
screen Eh . 

In short, the gap has a hybrid structure: The lower 
latitudes (with small h) are nearly vacuum having sub- 
GJ current densities and the inner boundary is located 
slightly inside of the null surface, because \& > region 
will be filled with the electrons emitted from the stellar 
surface. The higher latitudes, on the other hand, are non- 
vacuum having super-GJ current densities, and the inner 
boundary is located at the stellar surface, extracting ions 
at the rate such that their non-relativistic column den- 
sity at the stellar surface cancels the strong E\\ induced 
by the negative p— Pgj of relativistic electrons, positrons, 
and ions. The created pairs discharge such that E\\ virtu- 
ally vanishes in the region where pair creation is copious. 
Thus, in the intermediate latitudes between the sub-GJ 
and super-GJ regions, dp/ds > pa/ds holds. 

Even though the inner-most region of the gap is inac- 
tive, general relativistic effect (space-time dragging effect, 
in this case) is important to determine the ion emission 
rate from the stellar surface. For example, at h = h m /2 
for h m = 0.600 (i.e., the top right panel in fig. [3]), ji on is 
69% greater than what would be obtained in the Newto- 
nian limit, pgj = — £1- B /2ttc. This is because the reduced 
| Pgj | near the star (about 15% less than the Newtonian 
value) enhances the positive Eh , which has to be canceled 
by a greater ion emission (compared to the Newtonian 
value). The current, ji OI1 , is adjusted so that \p e s\ may 
balance with the trans-field derivative of \f near the star. 
The resultant |p e ff| becomes small compared to |pgj|j m 
the same manner as in traditional polar-cap models, which 
has a negative En with electron emission from the star. 
Although the non-relativistic ions have a large positive 
charge density very close to the star (within 10 cm from 
the surface), it cannot be resolved in figure [31 Note that 
the present calculation is performed from the stellar sur- 
face to the outer magnetosphere and does not contain a 
region with E\\ < 0. It follows that an accelerator having 
E\\ < (e.g., a polar-cap or a polar-slot-gap accelerator) 
cannot exist along the magnetic field lines that have an 
super-GJ current density created by the mechanism de- 
scribed in the present paper. 
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Fig. 3. Total (solid), created (dash-dotted), and 
Goldreich- Julian (dashed) charge densities in 
SlB{s,h)/(2irc) unit, for a; = 70° and h m = 0.060 
at nine transficld heights, h. If there is an ion emission 
from the stellar surface, the total charge density deviates 
from the created one. 
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